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I. INTRODUCTION 



Understanding the origin of the hierarchical structure of quark/lepton masses and their 
mixing angles is one of important issues in particle physics. Actually, several types of 
scenarios have been proposed so far. Among them, the democratic Ansatz on Yukawa 
matrices is an interesting approach in particular for the lepton sector, because 

it can lead to almost bi-maximal mixing angles. 

The exactly democratic form of the Yukawa matrix can be obtained by a certain flavor 
symmetry, e.g. 5*3 symmetry. However, the exact one by itself is not realistic, because it is 
a rank-one matrix. Thus, a certain pattern of small symmetry breaking terms are usually 
added by hand for the purpose to realize realistic Yukawa matrices. However, it seems that 
any comprehensive explanation leading to such pattern of symmetry breaking terms has not 
been given. 

Recently, two of the authors have considered dynamical realization of the almost demo- 
cratic Yukawa matrices p*, that is, each entry of the Yukawa matrix has the same infrared 
(IR) fixed point and the almost democratic Yukawa matrix can be realized by dynamics at 
the IR region whatever their initial values are. Actually, these fixed points correspond to the 
so-called Pendelton-Ross fixed point |7(. The exact IR fixed point is not realistic, because we 
just obtain the exact rank-one Yukawa matrix. Some region close to the IR fixed point would 
be interesting, and deviations from the fixed point correspond to small breaking parameters 
in the flavor symmetry approach to the democratic Ansatz. Thus, one of our purposes in 
this paper is to study more about renormalization group (RG) flows of Yukawa couplings 
in models with the above IR fixed points leading to the democratic Yukawa matrices. We 
shall show that we obtain certain patterns of RG-flows of Yukawa matrices, and study how 
much realistic Yukawa matrices for the quark and lepton sectors are derived from natural 
initial conditions, and which type of fine-tuning is required. 

Supersymmetric extension is attractive as new physics beyond the standard model. 
Within the framework of supersymmetric standard models, a flavor mechanism, which de- 
rives realistic quark/lepton masses and their mixing angles, would affect somehow their 
superpartners, that is, soft supersymmetry (SUSY) breaking squark/slepton masses and 
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trilinear couplings, the so-called A-terms. Although superpartners have not been detected 
yet, patterns of squark/slepton mass matrices and A-term matrices are strongly constrained 
by experiments on flavor changing neutral current (FCNC) processes That is, FCNCs 
require almost degenerate A-terms and soft scalar masses. In particular, the constraints due 
to the the K° — K° mixing and // — > decay are severe. In addition, CP- violating phases of 
A-terms are also strongly constrained by experiments on electric dipole moments (EDMs). 
These are the SUSY flavor and CP problems. A solution for theses problems is to realize 
the universal soft scalar masses and the universal A-terms, whose phases are aligned with 
the phase of gaugino masses, e.g. the gluino mass. In this case, the A-term has no physical 
CP-violating phase. Thus, it is quite important to examine these constraints on any SUSY 
flavor mechanism. 

So far various flavor symmetries have been also applied to suppress the flavor non- 
universality in SUSY breaking parameters, while the mass matrices are also explained by 
the same symmetries. As for the democratic mass matrices, the S3 symmetry also constrains 
structure of the SUSY breaking parameters However, the FCNC or the CP problem is 
not solved by the symmetry alone. Besides, the breaking effects of the flavor symmetry on 
the SUSY breaking parameters are unclear. 

It is known that when a model with one flavor has a Pendelton-Ross fixed point, the 
corresponding A-term A as well as the sum of soft scalar masses squared also has an IR 
fixed point like A = —M, where M is the gaugino mass Q] . It is expected that our 
model with three flavors has a similar fixed point, and A-terms are aligned dynamically at 
the IR region Indeed, we shall show IR fixed points of A-term matrices are universal 
Aij = —M. That is a favorable aspect, but the exact IR fixed point is not realistic, because 
we just obtain the rank-one Yukawa matrix on the exact IR fixed point as said above. Some 
region close to the fixed point would be interesting. Hence, in this paper we study RG-flows 
of Yukawa couplings and A-terms. We evaluate how much A-term matrices differ from the 
universal form, when Yukawa matrices are somehow realistic. 

This paper is organized as follows. In section II, we review briefly the Pendelton-Ross 
fixed point of Yukawa coupling and the corresponding A-term in a SUSY one-flavor model. 
In section III, we review briefly our three-flavor model, where all entries of a Yukawa matrix 
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3 



have fixed points. We study fixed points of the corresponding A-terms, Aij. In section IV, 
we study RG-flow of Yukawa couplings such as to obtain a certain pattern of the Yukawa 
matrix by natural initial conditions. We apply our results to the charge lepton sector. We 
also analyze RG-flows of A-terms and examine their non-universality. In section V, we 
study RG-flows of up and down sector Yukawa matrices. We use SU(5)' x 577(5)" model as 
a concrete model. We show that a certain pattern of up and down sector Yukawa matrices 
are obtained. We also discuss RG-flows of A-terms and consider FCNC and CP-violation 
constraints. In section VI, we comment on degeneracy of sfermion masses. Section VII 
is devoted to conclusion and discussion. In Appendix, we give examples with additional 
couplings, where values of fixed points are shifted significantly. 



II. IR FIXED POINT IN ONE-FLAVOR MODEL 
A. IR fixed point of Yukawa coupling 

First, we briefly review on the Pendelton-Ross IR fixed point in the model with a simple 
gauge group and a single Yukawa coupling, which corresponds to the superpotential, 

W = y$i$ 2 $3, (1) 

where $7s are chiral superfields. The one-loop RG equations for the gauge coupling g and 
the Yukawa coupling y are obtained as 

/i-j-^ = (aa y — ca g )a y , (3) 

where a g = g 2 / (8tt 2 ), a y = y 2 / (8ir 2 ), b is the one-loop beta function coefficient of the gauge 
coupling, and a and c are group-theoretical constants and both are always positive. Now, 
let us consider the RG equation of the ratio x = a y /a g , 

dx 

fi— = [ax — (c — b)]a g x. (4) 

iXiLL 

This equation has a nontrivial fixed point at 



The condition c — b > should be satisfied such that this fixed point is physical, i.e. x > 0. 

Actually, we can solve this equation exactly, and show this fixed point is IR attractive. 
However, analysis on linear perturbation around the fixed point x = x* + Ax would be 
useful for discussions in the following sections. Here we study the RG equation for Ax, 
which satisfies 

fi— — = (c — b)a g Ax. (6) 
Thus, this fixed point is IR attractive when c — b > 0. Actually, its solution is obtained as 

(b-c)/b 



Ax(fj) _ ( a g (n) \ c _ ( a g (/x) ' 
Ax(A) [a g (A)J ~U(A), 



(7) 



Large values of a g and (c — b) lead to stronger convergence toward the IR fixed point. 
Asymptotically non-free theories, i.e. b < 0, are favorable to realize strong convergence. 



B. Infrared fixed point of A-term 

Next, we review on fixed points of A-terms. In general, softly broken SUSY models have 
trilinear couplings of scalar components, that is, the so-called A-terms. In the above model, 
we would have the soft SUSY breaking term, 

^010203, (8) 

where (pi are scalar components of chiral superfields $j. Here and hereafter, we use the 
notation of A-term A = h/y. 

Now let us consider the RG equations of the gaugino mass M and the A-term. It is 
straightforward to obtain those RG equations in a simple model. However, it is convenient 
to use the spurion formalism to obtain RG equations of soft SUSY breaking terms in com- 
plicated models like models, which we shall discuss in the following sections. Such procedure 
is as follows. We replace a g and a y as [bj 

a g -> a g = a g (l + M6 2 + MO 2 ) + ■■■, (9) 
a y -> a y = a y {l - A9 2 - AO 2 ) + ■■■, (10) 

where is the Grassmann coordinate of the superspace and the ellipsis denotes terms with 
6 2 6 2 , which are irrelevant to the RG of the A-term, but relevant to soft scalar masses. The 
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couplings, dig and a y , satisfy the same RG equations as those of a g and a y . That is, the RG 
equations of the gaugino mass M and the A-term A corresponding to eqs.(J2J) and (j3J) are 
written as 

dM 

fi— = -ba g M, (11) 

(J/Uj 

dA 

/x— = aa y A + cMa g , (12) 

where a, b and c are the same constants as eqs.(j2J) and (J3J). We can find that AA = (A + M) 
has a fixed point in these equation. Actually, it satisfies the RG equation, 

fx— — = —aa g MAx + ax*a g AA. (13) 

Cl Li 



At x — x*, it reduces to 



= ax * ag ^A = (c - b)a g AA. (14) 



Note that the constant a is always positive. Thus, when the Yukawa coupling has the 
Pendelton-Ross IR fixed point, the corresponding A-term always has the fixed point, 

A = -M. (15) 

The deviation from the fixed point AA decreases in the same way as Ax satisfying eq. (6). 
Note that the gaugino mass M and the A-term A are complex parameters. Thus, the above 
relation, A = —M, is realized including their CP phases, that is, the CP phase of the A-term 
is aligned with one of the gaugino mass at the IR fixed point. 

III. INFRARED FIXED POINTS IN THREE FLAVOR MODEL 
A. Democratic fixed point 

In this section, we briefly review on the model with three flavors in which the demo- 
cratic form of the Yukawa matrix is realized dynamically by IR fixed point. 

We consider the SUSY model including three flavors of matter superfields Fi and /; 
{i = 1, 2, 3) and nine Higgs fields Hij with their superpotential, 

W = ^F l f J H lJ . (16) 

hi 



Anomalous dimensions of these fields, 7f;,7/ 3 - and 7^, are written as 



If, = a F oi Vik - c F a g , 
k=i 
3 



k=l 



lH i:j = 3a H a yi] - c H a g , 



(17) 

(18) 
(19) 



where a Vij = \yij\ 2 / (8n 2 ) , o-fj^h are positive constants and cfj,h are obtained as cfj,h = 
2C2(Rfj,h) by the corresponding quadratic Casimir C2(Rfj,h)- The RG equations of a yij 
are given as 



da 



//- 



d\i 



(20) 



As the previous section, we define Xij = a yij /a g . Then, it is found that the RG equations 
of have the fixed point, 



c-b 



XlJ ~ X ~ 3a ' 



(21) 



where a = ap + a/ + a H . To see whether this fixed point is IR attractive, we consider the 
linear perturbation around the fixed point, x^ = x* + Ax^. The RG equations of Axij are 
obtained as 



/i—Axij = 3aH<y g x*Axij + apa g x* Ax ik + afa Q x* Ax kj . 

k k 



(22) 



Alternatively, we can write more explicitly 



/i 



d 



dfi 



( A E E\ ( A^ 



A 2 
\A 3 J 

where Aj = (Axu, Axi2, Axis) T and 



E A E 
E E A 



A 2 

VA 3 y 



(23) 



A 



( a! dp Op ^ 

dp a' dp 
\clf clf a! J 



E 



I a f N 

a f 
V a j ) 



(24) 



where a' = dp + af + 3a H . This (9 x 9) matrix has the following eigenvalues, 



3a H , 3a H , 3(a F + a H ), 3a H , 3a H , 3(a F + a H ), 3(a f + a H ), 3(af + a H ), 3a. (25) 
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All of them are positive and it is found that the fixed point 0C%j — OC IS IR attractive. Hence, 
this model dynamically realizes the democratic form of the Yukawa matrix, 

(\ 1 l\ 
111 
1 1 1 



a yij oc 



V 



(26) 



/ 



To obtain the democratic form of the fermion mass matrix, we have to assume all of vacuum 
expectation values (VEVs) (Hij) are the same. Alternatively, we assume mass terms among 
Hij such that they lead to a single light mode H = 1/3 YUj Hij aric ^ the others have masses at 
a high energy scale. Here we take the latter scenario, and such mass terms have been shown 
in Refs. j^, [y]. We denote their mass scale by M H , and at M H the RG-flows approaching 
toward their fixed points are terminated. 

B. Fixed points of A-terms 



Here, we consider fixed points of A-terms. Using the spurion technique, we can write the 
RG equations of A^, 

/i— Aij = a F ^2 AikCtyik + a>f AkjOtykj + "iauAijOiyij + cMa g . (27) 
d ^ k k 

It is straightforward to show that there is the fixed point A^ = —M. Around the fixed 
points, we write the RG equations of AAij = A^ + M, 

d 



fji—AAy = a g x*{a F J2k ^A ik + a f J2k &A kj + 3a H AA i:j } 
-a g M{a F Efc ^Xik + afT,k ^%kj + 3a H Axij}. 



(28) 



When Axij = 0, AA^ satisfy the same RG equations as Axy. Thus, it is found that when 
the point IR attractive, the corresponding A-terms have the IR attractive fixed 

point, 



.4 



-M 



(29) 



1 1 l\ 
1 1 1 
1 1 1 

This form of A-terms is quite important from the phenomenological viewpoint. Obviously, 
the A-terms are universal, that is, favorable from the FCNC constraints. Furthermore, the 



V 



/ 
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CP phase of the A-terms is aligned with the CP phase of the gaugino mass M, that is, 
there is no physical CP-violation phase in the A-terms. Therefore, this form of A-terms 
can avoid both SUSY flavor and CP problems. Constraints from FCNC processes as well as 
CP- violations are usually presented by using the mass insertion parameters, e.g. (81 R )ij, 

„.a Aa 

= ~~2 ~ r " ■ (30) 

m SUSY 

in the S-CKM basis, where a = u,d,£ and msusY denote the average of sfermion masses 
and v a denotes the VEV of the corresponding Higgs field. 

IV. REALISTIC FERMION MASSES 
A. RG flow of Yukawa matrix 

In the previous section, we show the democratic form of the Yukawa matrix is realized at 
the IR fixed point, and at the same time the universal A-terms are obtained. That is quite 
favorable from the viewpoint of SUSY flavor and CP problems. However, the Yukawa matrix 
at the exact fixed point is not realistic, because it is a rank-one matrix, and only one family 
becomes massive and the other two remain massless. A region deviated slightly from the 
fixed point would be interesting. Thus, in this section we study the RG flows of the Yukawa 
matrix to investigate which pattern of the Yukawa matrix is obtained through a finite range 
of running by generic initial conditions and examine the possibility for realistic Yukawa 
matrices. In particular, we apply our analysis to the charged lepton sector. In addition, in 
a such situation, the A-terms may also deviate from the universal form somehow. We study 
whether such deviation leads to dangerous FCNCs. 

The democratic matrix can be diagonalized as follows, 



(l 


1 


l) 




f° 





0^ 






1 


1 


1 


= u 













(31) 


V 1 


1 


1 ) 




1° 












where 





( 


1 

V2 


1 




u = 




-1 
y/2 


1 

VE 


l 




\ 





-2 


vf / 



(32) 



9 



For example, if the lepton mixing matrix is obtained as Vmns ~ U T , that leads to almost 
realistic mixing angles, which are consistent with the neutrino oscillations 



3, i.e., 



MNS) 12 



0.5, 



MNS )23 



0.7, 



MNS ) 13 



< 0.1. 



(33) 



That is the reason why the democratic form of the Yukawa matrix is interesting, in par- 
ticular in the lepton sector. That is, the Yukawa matrix of charged leptons may almost be 
democratic, and it determines the mixing angles dominantly while the contribution from the 
neutrino Yukawa matrix to the mixing angle may be small 0. 

Sometimes this "hierarchical" basis is more convenient to discuss than the democratic 
basis. We denote the deviations of Xij from the fixed point in the hierarchical basis by Axij. 
Their RG equations are obtained as 



fj,—Axij(fj,) = a g x*Xi j Ax ij (fj l ), 



(34) 



with 



f 3a# 3clh 
y 3(a u + Off) 3(a u + a H ) 



Thus, they are solved as 



Axij(n) = Axij(A) 



+ clh) 
+ a H ) 
3a 



x \ 
b a *j 



(35) 



(36) 



and the Yukawa matrix in the hierarchical basis, yij is obtained as 



/ Axn (A) 



~ 3y/x* g(fi) 



6x 



3a H 



6.x 



6x* 
Ax 23 (A), 



£( n) 3 ( a Q+ a H) 
6x* ^ VA-V Qx* ^\t^J Qx* vrV 



(37) 



where the suppression factor e(n) is obtained as 

r a g (n)' 



-x*/b 



e(li) 



a g (A) j 



(38) 



Note that the powers of the suppression factor in the (i, 3) and (3, z) entries are larger 
than those in the entries = 1,2). Hence, when initial values Axij(A) are of the 
same order, e.g. of 0(1), the (i, 3) and (3, i) entries are suppressed rapidly and we obtain 



10 



the Yukawa matrix, 

{ l) 



jjij ~ 3Vx*5f(/i) 



(39) 



Thus, this RG-flow can explain why the third family is the heaviest. 

To obtain numerically realistic result of the lepton sector, we need e(fi) 3aa ~ m^/m T 
when Axij/(6x*) ~ 1. Furthermore, the (2,3) mixing angle is almost determined by the 
(2,3) entry of U T , because are suppressed rapidly. Thus, we obtain 

(Vmns)23 = (40) 

up to a contribution from the neutrino sector. This prediction is good when the contribution 
from the neutrino sector is not large 0. 

On the other hand, the above form can not explain the mass hierarchy between the first 
and second families when all of initial values Axij(A) are of the same order, e.g. of 0(1). If 
the initial values of the first and second families are chosen, e.g., 

Axu(A) Ax 12 (A)\ I , 

1 ' 1 ^ ] =Ax(A) 22 | _ V-M | ; (41) 

Ax 21 (A) Ax 22 (A) 

we obtain the correct mass ratio m e /m^. In the democratic basis, these initial couplings are 
given as 

' Ax n (A) Ax 12 (A)\ / 1 ~xfW 1 

= Ax(A) 22 _ (42 

Ax 21 (A) Ax 22 (A) ) \ 1 1 + ^5 

In addition, it also contributes to the mixing angle (Vmns)i2 as 



(V M ns)i2 ~ -^(1 - yrrie/Sm^), (43) 

up to the contribution from the neutrino sector. Indeed, the value (£7 T )i 2 = l/v^, which is 
derived from the exact democratic form, is slightly different from the value consistent with 
the neutrino oscillations. Thus, the above small modification is rather preferable. 
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B. RG flow of A-terms 



Here, let us study the RG-flow of the corresponding A-terms. In the hierarchical basis 
the RG equations (J28J) are written as 

fx — - — — = a g x* XijAAij — Ma g \ijAxij. (44) 

(Jj LL 

Their solutions are obtained as 

"Aiy(A) XijM(A) ( a a {jn)' 



Axij(fj). (45) 



A^-(A) b V a g( A ), 

Note that the damping factor is the same except the term including Axjj(/i)a 9 (/z)/a 9 (A). 

When we apply to the charged lepton sector, it is first noted that AAj3 and AA& for 
i,j = 1,2 are strongly suppressed as given in ([59]). Therefore, the branching ratios for 
the lepton flavor violating processes, r — > ^ and r — > ej, turn out to be much smaller than 
the experimental bounds. 

It is found that another lepton flavor violating process, \i — > ej, is also suppressed to 
some extent. Since it is expected that e(fi) 3aH ~ (m M /m r ), we may evaluate 

AAijifx) ~ (m /i /m T )F ii (A), (46) 
F y (A) = A Ay (A) - ^^A^.(A), (47) 

for z, j = 1, 2. Hence, the mass insertion parameters (^lr) 101 hj = 1> 2 are estimated as 

( 5 lk).. ~ (m,j,/m T )(m T /msusY), (48) 

where msusY is the average mass of sleptons, and we have taken i^j(A) = msusY- For 
example, in the case with tususy = 100 GeV, we obtain (S e LR )i2 = 10~ 4 . The second term in 
Fij(A) has more suppression factor like O(10 _1 — 10~ 2 ) when we take Ax 12 (A) = ^m e /m^ 
like eq. (}4~T|) . Furthermore, if \b\ = O(10), it would leads further suppression factor by 
(^(lO" 1 ). On the other hand, experimental bound on the /i — > e r y decay requires (S e LR )i2 < 
10 -6 . Furthermore, the experiment on the EDM of the electron requires Im(5 e LR )u < 
10~ 7 . If there is no suppression factor for the first term AAjj(A) in the -Fjj(A), then these 
constraints are not satisfied. Conversely, if we have another mechanism to suppress AAy(A) 
then the flavor non- universality generated through Ax is sufficiently small. 



* For example, a certain class of (string-motivated) supergravity models leads to Ay (A) = — M (A) as initial 
conditions. See e.g. Ref. [16| and references therein. 
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V. REALISTIC QUARK MASSES 



A. RG flow of Yukawa couplings 

In the previous section, we studied a single sector of a Yukawa matrix, and apply to the 
charged lepton sector. In this section, we study quark masses. For quarks, both the up and 
down sectors of Yukawa coupling matrices should have democratic fixed points. Otherwise, 
we can not realize small mixing angles. Thus, we have to extend the previous analysis to 
the case including the up and down sectors of Yukawa matrices. 

Here, we consider the concrete model, which has been proposed in Ref. |5j. It is the 
5*7(5)' x 5*7(5)" GUT model. We assume three families of quarks as well as leptons cor- 
respond to 10j and 5j for 5*7(5)' as the usual 5*7(5) GUT, and they are singlets under 
577(5)". We denote the gauge couplings g' and g" for 577(5)' and 5*7(5)", respectively. We 
assume that g' is strong, but g" is weak. We also assume that 5*7(5)' x 5*7(5)" is broken 
into the diagonal group 5*7(5) at/above the GUT scale. Then, we would obtain the usual 
5*7(5) GUT. The gauge coupling g of the diagonal group is obtained as 1/g 2 = l/g' 2 + l/g" 2 , 
and it is weak. The Pendelton-Ross IR fixed points can be realized for 5*7 (5)' within short 
running, because the gauge coupling g' is strong. That is one of the reasons why we extend 
the usual 577(5) GUT to the 577(5)' x 5*7(5)" GUT. § We assume that the mass scale M H 
is around the GUT scale. 

Another reason to consider the 5*7(5)' x 5*7(5)" GUT is related with the SUSY breaking 

parameters. In the asymptotically non-free gauge theories considered in the previous section, 

the running gaugino mass also decreases towards lower energy scale. Meanwhile the soft 

scalar masses are enhanced through the strong gauge interaction and become much larger 

than the gaugino mass. In the extended GUT, however, there are two gaugino masses M' 

and M" for the distinct gauge sectors of 5*7(5)' and 5*7(5)" respectively, and the gaugino 

mass M obtained after the symmetry breaking is given by 

M M' M" . , 

— = — + —• 49 

g 2 g,2 g,,2 

Since we suppose g' to be large, M is almost the same as M", which is not reduced from it's 



In extra dimensional models convergence toward IR fixed points may be sufficiently rapid [l^ in a single 
SU (5) model. Besides, the A-terms as well as the soft scalar masses are strongly aligned into the universal 
for asymptotically free gauge theories in the extra dimensions [lsj . 
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initial value due to smallness of the gauge coupling g". Therefore, the soft scalar masses, 
which may be enhanced to be of the order of the initial M", do not dominate over the 
gaugino mass M, as long as the gaugino masses M' and M" are given with their initial 
values of the same order. 

We introduce nine pairs of Higgs fields H u (5)ij and H d (5)ij, corresponding to the up and 
down sector Higgs fields. We have the following superpotential; 

W = y^lOilOjH^ij + //: / ; 10,5 ; //' , i5), ; , (50) 

The anomalous dimensions of matter fields and Higgs fields are obtained as 

3 

7io, = EK«*fc + a vki) + a d Q®tik} ~ a u Q Uyii - c Q a' g , (51) 
fe=i 

3 



75, 


k=i 

= 3a u H a u yij - c u H a' g , 


(52) 


1m. 

»j 


(53) 


»j 


= SajjOtyij — c d H a g , 


(54) 


6, 3ajj 


= 4, a u Q = 3, a d Q = 2, 


a d = 4, c Q = 36/5, c d = c H = 24/5, 



with 3a^ 

where a^ d = lutfl 2 / (Sir 2 ) and a' g = g' 2 /(8ir 2 ). Then, the RG equations of 0Cy..,a d .. are 
obtained as 

^ = (7io l +7io J +7^)<, (55) 
da d 

V-r 1 = (7io, + 75, + Inf) <■ (56) 



dfx 



ij 



We consider the RG equations of = ay'^/a' g , that is, 

+ E 4(4 + 4) - a Q x l - <>4K4-> ( 57 ) 
fc=i 



= {(b'-c Q -c d - 4) + 3a d H x 



3 



+ E + **) + E 4(4) + E - «K4> (58) 

fc=l fc=l k=l 

where b' is the one-loop beta function coefficient of g'. It is straightforward to show these 
RG equations have the following fixed points, 

X(a d + a d H + a d ) - 2a d a Y , 

x u _ x u* _ V Q < H < ) Q /p-gx 

tJ 10a%« + a d ) + 3a u H (a d Q + a d H + a d ) ' v ; 
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X 



VI 



X 



d* 



5a u Q X - (3a u H + Wa%)Y 



30a d Q a u Q - 3(3a u H + 10a%)« + a d Q + a d ) ' 



(60) 



where X = 2cq + c% — b' and Y = c Q + c d + c d H — b'. We consider perturbations around these 



fixed points, A 



A 6 - = A" + A" , 



It is convenient to define 



A, 



A" - A M -, 



to solve their RG equations. Actually, the RG equations of A^' ' are written as 
d 



(61) 



d 



Zalx u *a> g /±t p 
3 



/^ A ? = {E( 2 «Q( A ^ + A &) + H( A ^ + A i)) 



fc=i 



+3a« Ag - a£Ag - a u Q A^}a f g x u * , 

{E«A* + « + a d A d kj ) + 3<Ag - ^A*}a>*. 



(62) 

(63) 
(64) 



fc=i 



Easily, we can solve the RG equations of A^ as 



AJM = A<S(A) gg; 



(65) 



The other elements are mixed in the above equations. Here, we use the hierarchical basis. 
Through longsome but simple algebraic calculations, we find that the RG equations of Ag, 
Ag and A d ri for i,j = 1,2 are decoupled each other. Then we can solve them, 

-3a u H x u */b' 



AgO.) = Ag(A) (gg 



MM) = A J( A ) 

A 3^) = A 3,(A) 



-3a d H x d */b' 



-3(a d +a d H )x d * /b' 



(66) 
(67) 
(68) 



for i,j = 1, 2. However, the other elements A^ and Af 3 for i — 1,2,3 are mixed as 



b'a' 



da' 



^13 
^13 



x u *{3a u H + Aa u Q ) 



\ 



2x u *a u Q 



6a,QX d * 
3x d *(a d Q + a d H ) 



/A f 3 



V2 



A 



12? 



(69) 
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b'a 



d 



A s 

^23 
^23 



x u *(3a u H + Aa u Q ) 6a d Q x d * 
2x u *a u Q 3x d *(a d Q + a d H ) 



2V2 



(70) 



it , d 



( AS 
^33 

\ ^33 



x u *{3a u H + \Qa u Q ) 12a d Q x d * 

3x d *(a d Q + A d + a d H ) 

4 X 



2 X U Q 



UqX 



/ AS 

^33 

A d 
^33 



22 1 



(71) 



The solutions Af 3 (/i) and A^/x) are obtained as 



Af 3 (^/Af 3 (A) 
A d M/A d 3 (A) 



( 



R2 Q2 



( 



-n+/b' \ 



[R 2 r+ Q 2 r_) \(a g (»)/a> g (A)))-"-/ b ' 



+ 



KMK(A))- 3a ^ 76 ' 



key 



where n± are obtained as 



n± 



D + A± sJ(D - A) 2 + ABC 



with 



V 



A B \ I x u *{3a u H + Aa u Q ) Qa d Q x d * 
CD) { 2x u *a u Q 3x d *(a d Q + a d H ) 
and similarly r± are obtained as 



r± = 



D — A± J(D — A) 2 + ABC 



(72) 



(73) 



(74) 



(75) 



2B{-V) 

In addition, Ri and Qi are integral constants determined by initial conditions. Furthermore, 
the constants k\ and k\ are determined by solving the following equation; 



/ 1 _ A -B W ^ 



-C l-D 



h 2 
\k 2 j 



2 v / 2(-&') 



(Mi - a 



22;? 



(76) 



with I = —3afjX u */b'. Substituting explicit values, we obtain 

n± = 9x u * + 5x d * ± a/(9x m * - 5x d *) 2 + 72x u *x d *. 



(77) 
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It is found that both n± are positive, and obviously n + > n_. The solutions Af 3 (/i) and 
A^ 3 (/i) are the same as Af 3 (/i) and A^/i) except by replacing 



(A* - Af 2 ) 



A 5 



(78) 



Similarly we can solve Af 3 (//) and A 33 (/i), and they have rapidly damping behavior. Since 
their explicit forms are irrelevant to our discussions, we omit them. 

The solutions A^(/i) and Af 3 (//) for % = 1,2 include three damping factors. The most 
slowly damping term is important. Since we always have n_ < 6(= 3a#), we find that A^(/x) 
and A^ 3 (/i) behave like (a^)~ n -/ 6 '. Thus, the up and down-sectors of Yukawa matrices behave 
totally as 

e'(/i) n - 




/ e '(^)3ai*<* 





1 



(79) 



£'(/^) n - 
1 



(80) 



where 



-1/6' 



(81) 



We have omitted the coefficients. Also we have omitted (3,1) and (3,2) entries in the down- 
sector Yukawa matrix, because they are damping more rapidly than the other entries in 
the down-sector Yukawa matrix. That is the same behavior as in the case with a single 
Yukawa matrix, which discussed in the previous section. However, (1, 3) and (2, 3) entries 
are damping not rapidly in the case including both the up and down-sector Yukawa matrices. 
That would be important to derive the mixing angle V c b in the quark sector. 

Recall that 3a^ = 6 and 3a^ = 4. Explicit values of n_ are n_ = 3.5,4.6 and 5.3 for 
x u * = 0.5x d *,x d * and 2x d *, respectively. Thus, when x u * ~ x d * , we can obtain the realistic 
relations, 



m t \m b J 

However, the model which we are discussing leads to 



Vcb ~ • 



x" 



255 306 



x' 



d* 



44 7 
— - —b', 
85 204 ' 



(82) 



(83) 
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and it predicts rather large ratio of m c /m t like m c /m t ~ ra s /m&. However, fixed point 
values x u * and x d * change significantly when we include other couplings. In Appendix, we 
give examples, where fixed points are shifted by adding additional couplings. In a such case, 
the above analysis is the same except by using new fixed point values x u * and x d * when 
additional couplings are close to their fixed points. Thus, we could obtain models with 
several couplings leading to the realistic mass relations of the second and third families, e.g. 
the model leading to x u * ~ x d *. Here we do not study explicitly a concrete model with 
values of x u * and x d * leading to realistic results. We use x u * and x d * as free parameters to 
present generic models with several couplings and effects due to such additional couplings. 

We have shown that we can realize the quark mass ratios and the mixing angles between 
the second and third families when fixed point values x u * and x d * are in the proper region. 
However, only the RG dynamics can not lead to the realistic mass hierarchy between the 
first and second families. That is the same as the situation in the previous section. Thus, 
we need fine-tuning of initial conditions. For example, if the initial conditions are chosen, 
e.g., 



( A u 
^11 


A" 

^12 


A" 

\ ^21 


A" 

^22 




A d 

^12 


K Z4i 


A d 

^22 



A M (A) 







\ 



22 



A d (A) 22 



(84) 



(85) 



/ 



we obtain the correct mass ratios m u /m c and md/m s , and the mixing angle V us ~ m^/m s . 
These initial conditions are also consistent with those in the lepton sector, which was dis- 
cussed in the previous section, because m e /m^ ~ m^/m s . 

B. A-terms 



Here, we study the corresponding A-terms. By using the spurion technique, we can obtain 
the RG equations of A^f and find they have the fixed points, 



AU,d 



A* = -M. 



(86) 



We expand A?f around the fixed point as A™f = A* + AA?f, and the deviations A A"/ 
satisfy the following RG equations; 



„—AA s 



a/ [3a^"*AAg - a^*AA* - a u Q x u *AA^ 
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+ £{2a^*(A^ + AA s kj ) + 2a d Q x d *(AA d k + AA d k )}} 
k=i 

-ajM[3a£Ag - a£Ag - a^Af 



+ £{2a£(Af fc + A* ) + 2a d Q (A d k + A d k )}], (87) 



fc=i 



/^ A 4 = «;[3a^x rf *A4 - ^x w A^g 



+ E{^"* A 4 + «^*A^ + a d x d *AA d kj ) }] 
k=i 

-ajM{3a£Aj. - |> A* + £(a&A* + ajA& + A'AjJ}, 



fc=i 



^|-A4 = 3alx u *AA? ja > g - a' g M3a u H A? p (89) 

where AAf- = AAf- + A A u -- and AA£ = A Af- - AA%. Thus, we see that AA^ behave like 
AAij ~ Aij similarly to the model in the previous section. 

The experimental bounds for AMk, AMb and the branching ratio of b — > S7 restrict the 
mass insertion parameters as (5^) 12 < 4.4 x 1CT 3 , (5£r) < 3.3 x 1CT 2 and (Sf^j < 
1.6 x 10~ 2 respectively. On the other hand, we can estimate the mass insertion parameters 
of the present model (S^n) for i,j = 1, 2 as 

V / ij 



(^LRji - ~ ( m s/msusY), (90) 



where we have taken AAy(A) = msusY- Then we obtain = 1CT 4 , which is consistent 

with the experiment on AMk- We have taken msusY = 500GeV. The mass insertion 
parameters (5f^j for i,j = 1,2 are also found to be comparable with (5lr) 12 , since 
the parameter n_ introduced in the previous subsection is not so different from ?>a d H x d * . 
Therefore, the bounds from AM B and b — > S7 are also satisfied. 

Moreover, the EDM of the neutron requires Im ((5l'r) ) = 1CT 6 , which seems to be 
rather severe in general. However, this condition may be also explained, when \Aij\ ~ M' is 
smaller than M" ~ M by one or two orders at the GUT scale. This is because the A-terms 
become almost flavor universal at the weak scale due to corrections through the MSSM 
gauge interactions, which dominate over the initial values. Otherwise, we need fine-tuning 
of initial conditions AAij (A) or some mechanism to lead to such required initial conditions. 
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VI. COMMENT ON SFERMION MASSES 



In the previous sections, we have concentrated to constraints on the A-terms. Indeed, we 
have stronger constraints on the A-terms from FCNC and CP processes than soft sfermion 
masses. Here we give a comment on sfermion masses. From the viewpoint of experiments 
on FCNC processes, degenerate sfermion masses are favorable. 

The dynamics due to IR fixed points leads to the certain relation among soft scalar 
masses. When the Yukawa coupling yijFifjHij has an IR fixed point, the sum of soft scalar 
masses squared also has the IR fixed point, 

m\ + m).+m 2 Hi . = \M\ 2 . (91) 

However, these relations are not enough to lead to degenerate sfermion masses among flavors 
when there are nine Higgs fields H^. If there is any symmetry relating the soft scalar masses 
as well as the Yukawa couplings of the nine Higgs fields, then the above relations would lead 
to degenerate sfermion masses. An example is the model with symmetry 19j . In this 
model, three families correspond to triplets under A4, and four pairs of Higgs fields are 
introduced, and they correspond to A 4 singlet and triplet. That is, three Ha are identified, 
i.e., Hn = H 2 2 = H 33 , and ify for i 7^ j is identified with Hji, = Hji, in words of 
our model with nine Higgs fields. In such model, degeneracy of sfermion masses can be 
realized as IR fixed points [20j. However, in this model entries of the Yukawa matrix are 
also identified by the A^ symmetry, that is, all of diagonal entries yu are always the same, 
i.e., yn = y 2 2 = 2/33, and all of off-diagonal entries y^ (i 7^ j) are always the same, i.e., 
Vij = Vjii even away IR fixed points. We would obtain realistic results from such model 
when yu = y^ for i 7^ j and VEVs of four pairs of Higgs fields are fine-tuned in a proper 



way 
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21| . Because that is a scenario different from the scenario in this paper, we do not 



study it here, but it may be another interesting scenario. 



VII. CONCLUSION AND DISCUSSION 



We have considered the model to realize the democratic form of Yukawa matrices. We 
have studied the RG flows of Yukawa couplings and obtained their specific patterns. The 
mass hierarchy and the mixing angles between the second and third families can be realized. 
For a single sector of the Yukawa matrix, the mixing angle between the second and third 
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families is determined almost by U% 3 . However, for the up and down sectors of Yukawa 
matrices the mixing angle between the second and third families can be obtained as a value 
similar to m s /mf,, i.e., m s /m;, ~ V^. These aspects are quite interesting. Also, we can 
obtain the mass hierarchy, m s jm^ ~ m^/m T and m c /m t ~ (m s /mb) 3 ^ 2 . 

The corresponding A-term couplings also have the universal IR fixed point. That is also 
important in order to avoid SUSY FCNC and CP problems. In particular, there are strong 
constraints from the Kaon system, the fi — > c-f and EDMs. FCNC constraints can be relaxed 
in the parameter region realizing realistic values of quark mass hierarchies, except for the 
— > e7 process. Constraints due to CP violations can be ameliorated, but that may not be 
sufficient. We would need fine-tuning of the initial conditions on the A-terms. Alternatively, 
we may need some mechanism to realize such fine-tuning. 

Our model does not intend to explain the mass hierarchy and mixing angle between the 

first and second families. In order to do this, we need to choose proper initial conditions 

for them, or to find some mechanism to realize such initial conditions. For example, we 

may consider to achieve dynamical alignment of Yukawa couplings so that the large mass 

fl 

hierarchy like m u /m t can be generated [5J]. This can be realized by considering a strongly 
coupled GUT. Moreover, it was shown in Ref. |5j that the mass matrix close to the Fritzsch- 
type one can be obtained, if one of the Yukawa couplings is given with a relatively small 
initial value. In such a situation, the gauge dynamics aligns the A-terms strongly enough, 
and the problematic processes of FCNC and CP violations may be automatically suppressed 
well below their experimental bounds. That is beyond the scope of this paper and we would 
study it elsewhere. 

Within the framework of our model with nine pairs of Higgs fields, we can not control 
squark/slepton masses only by fixed point dynamics such that they are degenerate to avoid 
FCNC constraints. 
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Appendix 

Here we give two examples of models with additional couplings, which shift values of 
fixed points. 

First, we consider the model, which has the coupling among H^- d and the adjoint field 
24. The relevant superpotential is written as 

W = y^WWH u {5) l0 + / / ; , ; l() ; ry//' , (5!, ; + Kij, kl H u (5) ij H d (5) kl 24. (92) 

There may be other couplings as well as other fields. However, we assume that only the 
above couplings are around their fixed points, and we neglect other couplings. In this model, 
anomalous dimensions are obtained as 

3 

7i 0l = EK« fe + «L) + «J - «Q< - cg«i, (93) 
k=i 

78, = « d E<-Q<, (94) 
k=i 

7fl}$ = 3a H a y ZJ +d H J2 a vM - c u H a' g , (95) 

k,l 

l H d = ^ a H ad yij +d H J2 a kl,ij ~ C H a 'gi ( 96 ) 

13 k,l 

724 = a 24 E a nM - c 24<y' g , (97) 

ij,kl 

where a 2 4 = 1, c 24 = 10, dn — 5 and a^^i = |/tij,fc/| 2 /(8vr 2 ). The RG equations of Yukawa 
couplings are written as 

= (7io, + 7io 3 + 7^)<, (98) 



d 

/i— a 
dji 



■ d yij = (7i 0i +75,+7^)<, (99) 



V-^ a iJM = (7flg + 7n* + 724)aij,fci- (100) 



This model has the IR fixed points, 

3959 - 1756' 



x u * = ^[ 363 f = °- 29 " 0-0136', (101) 
* d * = 3361 90~9Q 245& = 0-37- 0.027^, (102) 
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Next we consider the model, which includes extra fields 5' and it couples with Hfj and 
the adjoint field. Other extra fields are added such that this model is anomaly-free. The 
relevant superpotential is written as 

W = ^10 l 10^ u (5) 4j + // ( ''/10'ry7/' , (5), ; + K' kl 5'H d (5) kl 24. (103) 

The anomalous dimensions are obtained as 

7io i; = EK« fe + <,) + «J - - c Q a' g} (104) 

k=l 

75, = a d E<-c d «;, (105) 

k=l 

lH u = 3a%al. - c u H a' g} (106) 

lHf . = 3a d H a d yij + dua'^ - c d H a' g (107) 

724 = diA^a'ij - c 2 4«3, (108) 

ij 

75' = d H J2 a 'ij - c Ha'g, (109) 

ij 

where a'^ = |Kjj| 2 /(87r 2 ). We can obtain the RG equations of Yukawa couplings. Then we 
find that they have the IR fixed points, 

12828 - 5756' 



x 



0.47- 0.0216', (110) 



27270 
3432 -365* 

18180 V ; 

Thus, as seen in these two models, values of fixed points x u * and x d * are shifted significantly. 



[1] H. Fritzsch, Nucl. Phys. B 155, 189 (1979); Y. Koide, Phys. Rev. D 28, 252 (1983); Phys. 
Rev. D 39, 1391 (1989). 

[2] H. Fritzsch and Z. Z. Xing, Phys. Lett. B 372, 265 (1996) |arXiv:hep-ph/ 9509389 ; 
Phys. Lett. B 440, 313 (1998) |arXiv:hep- ph/9808272 ; Phys. Rev. D 61, 073016 (2000) 
|arXiv:hep-ph/9909304| ; Phys. Lett. B 598, 237 (2004) |arXiv:hep-ph/0406206| . 

[3] M. Tanimoto, Phys. Lett. B 483, 417 (2000) |arXiv:hep-ph/0001306| ; M. Fukugita, M. Tani- 
moto and T. Yanagida, Phys. Rev. D 57, 4429 (1998) |arXiv:hep-ph7 9709388 . 

[4] M. Tanimoto, T. Watari and T. Yanagida, Phys. Lett. B 461, 345 (1999) 
arXiv:hep-ph/9904338 ; E. K. Akhmedov, G. C. Branco, F. R. Joaquim and J. I. Silva-Marcos, 

23 



Phys. Lett. B 498, 237 (2001) |arXiv:hep-p"h /0008010 ; G. C. Branco and J. I. Silva-Marcos, 
Phys. Lett. B 526, 104 (2002) |arXiv:hep-ph/0106125| ; T. Watari and T. Yanagida, Phys. 



Lett. B 544, 167 (2002) arXiv:hep-ph/0205090 ; Q. Shan and Z. Tavartkiladze, Phys. Lett. 



B 594, 177 (2004) ar Xiv:hep-ph/ 0401235 . 
[5] T. Kobayashi, H. Shirano and H. Terao, Prog. Thoer. Phys. 113 (2005) 1077 

|arXiv:hep-ph/0412299] . 
[6] S. A. Abel and S. F. King, Phys. Lett. B 435, 73 (1998) |arXiv:hep-ph/9804446| . 
[7] B. Pendleton and G. G. Ross, Phys. Lett. B 98, 291 (1981). 

[8] F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, B 477, 321 (1996) 
|arXiv:hep-ph/9604387] . 



[9] P. H. Chankowski, O. Lebedev and S. Pokorski, arXiv:hep-ph/0502076 
[10] K. Hamaguchi, M. Kakizaki and M. Yamaguchi, Phys. Rev. D68, 056007 (2003) 
|arXiv:hep-ph/0212T72] . 

[11] M. Lanzagorta and G. G. Ross, Phys. Lett. B 364, 163 (1995) |arXiv:hep-ph/9507366| ; 
P. M. Ferreira, I. Jack and D. R. T. Jones, Phys. Lett. B 357, 359 (1995) 
|arXiv:hep-ph/9506467] . 

[12] T. Kobayashi and K. Yoshioka, Phys. Lett. B 486, 223 (2000) |arXiv:hep-ph/0004175| ; 
T. Kobayashi and H. Terao, Phys. Lett. B 489, 233 (2000) |arXiv:hep-p"h7 0005265 . 

[13] T. Kobayashi and K. Yoshioka, Phys. Rev. D 62, 115003 (2000) |arXiv:hep-ph/0005009| . 

[14] Y. Yamada, Phys. Rev. D 50, 3537 (1994) |arXiv:hep-ph/9401241| ; J. Hisano and M. A. Shif- 
man, Phys. Rev. D 56, 5475 (1997) |arXiv:hep-ph/9705417| ; I. Jack and D. R. T. Jones, 
Phys. Lett. B 415, 383 (1997) |arXiv:hep-ph /9709364 ; I. Jack, D. R. T. Jones and A. Pick- 
ering, Phys. Lett. B 426, 73 (1998) |arXiv:hep-ph/9712542| ; Phys. Lett. B 432, 114 (1998) 
|arXiv:hep-ph/9803405] ; L. V. Avdeev, D. I. Kazakov and I. N. Kondrashuk, Nucl. Phys. B 
510, 289 (1998) |arXiv:hep-ph/9709397| ; T. Kobayashi, J. Kubo and G. Zoupanos, Phys. Lett. 



B 427, 291 (1998) arXiv:hep-ph/9802267 ; N. Arkani-Hamed, G. F. Giudice, M. A. Luty 
and R. Rattazzi, Phys. Rev. D 58, 115005 (1998) |arXiv:hep-p"h"79 803290 ; D. I. Kazakov 
and V. N. Velizhanin, Phys. Lett. B 485, 393 (2000) |arXiv:hep- ph/0005185 ; H. Terao, 
arXiv:hep-ph/0112021| 



[15] M. Maltoni, T. Schwetz, M. A. Tortola and J. W. F. Valle, New J. Phys. 6, 122 (2004) 



arXiv:hep-ph/0405172 . 



24 



[16] Y. Kawamura, T. Kobayashi and J. Kubo, Phys. Lett. B 405, 64 (1997) 
|arXiv:hep-ph/9703320] . 

[17] M. Bando, T. Kobayashi, T. Noguchi and K. Yoshioka, Phys. Lett. B 480, 187 (2000) 
|arXiv:hep-ph/0002102| ; Phys. Rev. D 63, 113017 (2001) |arXiv:hep-ph/0008120| . 

[18] J. Kubo and H. Terao, Phys. Rev. D 66, 116003 (2002); Y. Kajiyama, J. Kubo and H. Terao, 
Phys. Rev. D 69, 116006 (2004) |arXiv:hep-ph/0311316|| ; K-Y. Choi, Y. Kajiyama, H.M. Lee 
and J. Kubo, Phys. Rev. D 70, 055004 (2004) |arXiv:hep-ph/0402026| . 

[19] K. S. Babu, T. Enkhbat and I. Gogoladze, Phys. Lett. B 555, 238 (2003) 



arXiv:hep-ph/0204246 . 

[20] K. S. Babu, T. Kobayashi and J. Kubo, Phys. Rev. D 67, 075018 (2003) 



|arXiv:hep-ph /0212350 . 
[21] P. M. Fishbane and P. Kaus, Phys. Rev. D 49, 4780 (1994); Z. Phys. C 75 (1997) 1. 



25 



